Magnetotransport near a quantum critical point in a simple metal 
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We use geometric considerations to study transport properties, such as the conductivity and Hall 
coeflicient, near the onset of a nesting-driven spin density wave in a simple metaL In particular, 
motivated by recent experiments on vanadium-doped chromium, we study the variation of transport 
coefficients with the onset of magnetism within a mean-field treatment of a model that contains 
nearly nested electron and hole Fermi surfaces. We show that most transport coefficients display 
a leading dependence that is linear in the energy gap. The coefficient of the linear term, though, 
can be small. In particular, we find that the Hall conductivity axy is essentially unchanged, due to 
electron-hole compensation, as the system goes through the quantum critical point. This conclusion 
extends a similar observation we made earlier for the case of completely flat Fermi surfaces to the 
immediate vicinity of the quantum critical point where nesting is present but not perfect. 

PACS numbers; 72.15.Eb, 71.10.Hf, 71.18.+y, 75.10.Lp 
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I. INTRODUCTION 

Continuous zero temperature phase transitions have 
attracted much interest lately, due to a variety of un- 
usual low-temperature properties observed in their vicin- 
ity. Experiments on a broad range of materials, including 
cuprates, heavy fermion compounds and transition met- 
als, have shown thermodynamic, magnetic and transport 
properties that often do not fit in the framework of the 
Landau Fermi liquid theory of metals. This raises the 
possibility that the critical ground state may be a new, 
exotic state of matter. In many cases, understanding 
such a transition is rendered difficult by the fact that the 
phases involved are not well understood. This is the case 
for the cuprate superconductors, where even the under- 
lying 'normal' state remains a subject of some debate, 
and for many heavy fermion materials, which involve a 
nontrivial interplay between the f-electron magnetic mo- 
ments and conduction electron sea. Thus, it is important 
to study a case of a simple magnetic metal, where both 
the magnetic and the paramagnetic phases, as well as the 
nature of the transition, are believed to be well under- 
stood. 

The Cti-x'Vx alloy is an ideal candidate for such a 
study (see reviewai^). Pure chromium is a textbook ex- 
ample of a nesting-driven antiferromagnet, where doping 
with vanadium suppresses the Neel temperature T^r to 
zero already at a; = ~ 3.5%. In recent experiments 
on Cri_a;Va;^ Tat has been driven to zero by pressure, 
which allows access to the quantum critical point with 
minimal involvement of doping-induced disorder. 

Of various thermodynamic and transport properties 
near a quantum critical point, the Hall transport at- 
tracted particular interest lately, as it carries informa- 
tion about the change of the Fermi surface volume in an 
interacting system^. Here we will stay within an SDW 
description, for which the corresponding quantum critical 



point is above the upper critical dimension&i. We will 
focus exclusively on the static transport coefficients and, 
with the goal of treating the problem within the sim- 
plest possible model, confine ourselves to a mean-field 
approach. While the behavior of the electric conduc- 
tivity near the onset of a density wave within such an 
approach has been studied early on^, the Hall conductiv- 
ity has not been investigated until recently, when it was 
discussed for an SDWS and for d-density wave order—. 
Standard expressions^^ give electric and Hall conductiv- 
ities in terms of integrals over the Fermi surface. When 
evaluating them, one has to take into account the Fermi 
surface change in the ordered state, as well as the change 
of Fermi velocities and mass tensors. Consequently, de- 
veloping intuition about the behavior of these quantities 
is non-trivial. 

In our earlier work^, we considered a model in which 
the Fermi surface contains completely flat portions. In 
that case, the Hall number has a discontinuity at the 
quantum critical point, but the Hall conductivity does 
not change due to the zero curvature of the flat parts. 
Moreover, if the tuning parameter (such as pressure) 
does not change the elastic scattering, the Hall coefficient 
should scale as the square of the longitudinal resistivity. 
This behavior was confirmed by detailed numerical sim- 
ulations based on the actual Fermi surface of chromium, 
and the results were in good agreement with Hall data— 
across the doping induced quantum phase transition in 
Cri_,V,. 

Fermi surfaces, however, are never perfectly flat. The 
finite curvature becomes especially important in the im- 
mediate vicinity of the quantum critical point, which 
has become experimentally accessible with pressure in 
Cti-x^x^ In this paper, we address this regime in some 
detail by considering the Rice modeU^. The model, mo- 
tivated by the electronic structure of chromium, incorpo- 
rates an electron Fermi surface and a hole Fermi surface 
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that are nearly nested but do have a finite curvature; 
moreover, it takes into account an electron reservoir that 
mimics other bands. We determine the transport coef- 
ficients using a geometric approach. We find that most 
transport coefhcients display a leading dependence that is 
linear in the energy gap. However, due to cancellation of 
the electron and hole contributions, the Hall conductivity 
axy is essentially unaffected by the onset of magnetism, 
even in the immediate vicinity of the quantum critical 
point. The insensitivity of a^y to magnetism is, therefore, 
a robust feature of the model over a wide range of pa- 
rameters across the magnetic quantum phase transition. 
This conclusion is in agreement with the experiment a^i^^ . 

The remainder of the paper is organized as follows: 
In Section H, we discuss the Fermi surface of Cr, and 
the Rice model used in this papeiii. In Section HI, 
we discuss the changes in the electron spectrum upon 
magnetic ordering. In Section IV, we describe a geo- 
metric pictureii that affords a simple description of the 
transport coefficients near the quantum critical point in 
Cti-x^x- In Section V, we discuss numerical calculations 
based on the actual band structure of Cr. Some conclu- 
sions are offered in Section VI. Calculational details are 
described in the Appendices. 



II. MODEL FERMI SURFACE OF CHROMIUM 

The Fermi surface of chromium consists of several 
sheets and is well accounted for by band structure 
calculations^'^. It is believed that the interaction be- 
tween two particular sheets: the F-sheet, centered around 
the F-point, and the H-sheet, centered around the H- 
point, is responsible for the magnetism. These sheets are 
nearly nested; as a result, an SDW with nesting wave 
vector Q emerges, as shown in Fig. ^ In chromium, the 
SDW is of the linear type, which means^ that both +Q 
and — Q nesting is relevant. 

The relative size of the Fermi surface (FS) sheets can 
be changed by doping with vanadium. It is believed that 
the dominant effect of the doping is to change the elec- 
tron concentration, which corresponds to the shift of the 
chemical potential ^. As /i shifts with increasing hole 
doping X, the mismatch between the F-surface and the 
H-surface increases. The H-sheet grows and the F-sheet 
diminishes in size. This leads to the change of wave vec- 
tor with concentration Q — Q(a;). At the critical doping 
Xc, the mismatch becomes so large that the SDW forma- 
tion is not energetically favorable any more. 

To allow for an analytic treatment, we will employ 
the Rice model of itinerant antiferromagnetismi^ that 
involves an electron and a hole Fermi surface of spherical 
shape and unequal size. The Fermi surface of chromium 
is certainly far from being spherical^. However, it is not 
completely flat, and therefore the nearly nested parts can 
be fairly accurately approximated by spherical segments. 
It is these parts that are important in the development 
of the SDW. Also, as long as the SDW gap is small, only 
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FIG. 1: Top: a 2D cross section of the relevant Fermi sur- 
face sheets in Cr. The electron-like Fermi surface centered 
at the F-point and the hole-like Fermi surface centered at 
the H-point have regions that are almost nested. These re- 
gions lead to the SDW formation. The dashed lines show 
the H-centered Fermi surface displaced by the vectors ±Q. 
Bottom: schematic view of the electron and the hole energy 
bands along the F-H line. The dashed lines shows the hole 
band displaced by ±Q. 



these parts are destroyed by the gap opening. This is 
why a model that considers a spherical FS will give the 
correct behavior of physical quantities as A 0. The 
highly curved, non-nested parts of the F-and H-sheets 
will introduce corrections only when A becomes large. 

We consider spherical Fermi surfaces of radii kpe for 
the electron surface and kph for the hole surface. Wher- 
ever possible [for exceptions, see Sees. IIV A2llTVB 2| . we 
will linearize the spectrum: 



f-h{k) = -Vh{\k\ - kph) 



(1) 



We will present an analytic treatment of two cases (see 
Fig. 121). In the first case, upon translation by Q, the 
electron Fermi surface crosses the hole Fermi surface. In 
the second case, the two Fermi surfaces touch each other 
upon translation. The latter case is never realized in the 
Rice model: the maximum in the susceptibility is at a Q 
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FIG. 2: The electron-like FS (solid line) and the hole-like FS 
displaced by -|-Q (dashed line). The arrows show velocities 
at the Fermi surface. They point outside the electron Fermi 
surface and inside the hole Fermi surface. (A) the "crossing" 
case, (B) the "touching" case 



vector corresponding to crossing spherea^^. 

Note that in Fig. |21 the model hole surface is shown 
displaced only by the -l-Q vector, while in the real sit- 
uation of chromium (Fig. ^ both the -t-Q and the — Q 
displacements lead to nesting and should be considered. 
This is done for simplicity, and also because our analytic 
calculation is intended to give only the functional form 
of the A dependence. Taking into account both displace- 
ments changes only the value of the coefficient, not the 
leading power of the dependence. In Appendix^ we ex- 
plain how the full result for the linear SDW is related to 
that of the -l-Q displacement only. 

There are other sheets of the Fermi surface of 
chromium that are not shown in Fig. ^ They do not 
cause the SDW formation, but are also modified by the 
onset of the SDW, and thus their contribution to the 
transport coefficients changes as well. The changes of 
those Fermi surfaces will correspond to the crossing case 
of Fig. 13 and our analysis applies to them as well, as 
long as we are interested only in the leading power in the 
dependence of transport coefficients on the energy gap 
A. 



III. SPECTRUM IN THE ORDERED STATE 

Onset of a linear SDW generates a non-zero matrix ele- 
ment A between the single-electron states connected by 
the ordering wave vector ±Q. These matrix elements are 
independent of spin* (as opposed to the case of a helical 



SDW). As a result, a gap opens at all momenta k and 
k±Q, related by the resonance condition e(k) — e(k±Q), 
where e(k) is the single-particle spectrum in the param- 
agnetic state. When the chemical potential falls inside 
the gap, some parts of the original (paramagnetic state) 
Fermi surface disappear, and the FS sheets reconnect as 
shown in Fig. |21 

In the following, we neglect smaller gaps that appear 
whenever two single- particle states meet a higher-order 
resonance condition e(k) = e(k -I- nQ) with n > 1, 
since the magnitude of these higher-order gaps decreases 
rapidly with n, as A„ ~ A(A/ei?)"^^- 

In chromium, the condition e(k) — e(k± Q) is fulfilled 
only for electrons belonging to different (hole-like and 
electron-like) bands. The same is true for the Rice model. 
Thus the spectrum in the ordered state takes the form 



61 (k) = 

e2(k) = 




(2) 



As discussed in the previous section, to find the func- 
tional dependence of physical quantities on A, we will 
use only one sign of Q in the formula for t2 (see Ap- 
pendix^. Thus we will consider one electron band and 
one displaced hole band. 

The original spherical Fermi surfaces of electrons and 
holes are given by the equations 

ee,/^(k) = 0. 

If the electron reservoir provided by other sheets of the 
Fermi surface is large, the gap opening does not lead to 
a shift of the chemical potential. The new Fermi surfaces 
are then given by 



e±(k) = 0, 



(3) 



As shown in Fig. |2 when the SDW gap opens, the Fermi 
surfaces reconnect. Instead of the two original spherical 
Fermi surfaces, one gets a larger hole FS and a smaller 
lens-shaped electron FS. The latter corresponds to the 
{+) sign, while the former corresponds to the (— ) sign in 
Eq. 13 In the touching case, there is no lens, thus only 
the (— ) sign is relevant. For either sign, this equation is 
equivalent to 



ei(k)e2(k) = A2. 



(4) 



This may be complicated in fc-space, but is very simple 
in the space of (ei, £2). 

The velocities in the modified bands are: 



1 d€± _ 
h~dk ' 

Vi -I- V2 



(5) 



± 



ei - £2 



Vi - V2 



£2 



4A2 



4 



In this formula, velocities on the right hand side are 
taken at the same momentum k, so one or both of them 
have to be evaluated far from the corresponding original 
electron or hole Fermi surfaces. 

Away from the points where the electron and the hole 
Fermi surfaces cross or touch each other, the new Fermi 
surface is close to one of them. Quantitatively, this means 
that the distance between the new and old surfaces is 
much smaller then the distance between FSi and FS2. 
The corresponding condition reads (for definiteness, we 
write it here for one of the branches): 

|e+ - eil << |ei - £2], 

which, using can be re-expressed as 

|ei -£2! » A. 

The crossover from large to small deviation of the new 
Fermi surface from the old one takes place at 

lei-esi ~ A. (6) 

The above equations neglect the variation of the chemi- 
cal potential upon onset of ordering. The ordering wave 
vector Q also changes, and the dependences A(a;), /i(a;) 
and Q{x) on doping, x, need to be found self-consistently. 
Within the framework of the Rice model, this was done 
within a mean-field treatment in the original paper—. 
Both ^i{x) and Q{x) turn out to depend on doping lin- 
early: 

fi{x) = ^{xc) + ^J.' -{x-Xc), Q{x) = Q{xc) + Q' -{x-Xc), 

while A(a;) shows the standard mean-field dependence 
A(a;) ~ \/xc — X. Hence, in the immediate vicinity of 
the transition, the leading dependence of the transport 
coefficients on doping is defined implicitly by the depen- 
dence on A(x), and can be found neglecting the variation 
of /i and Q. We now turn to finding this dependence. 

IV. TRANSPORT COEFFICIENTS 

A. Conductivity 

In the Boltzmann equation approach, the conductivity 
is given by the integral 

eV f dSp 

where Va — {l/h){de±/dka), e is the electron charge, 
dSp is the Fermi surface element in fc-space, and r is the 
electron scattering time, hereafter assumed isotropic^'^. 
In the SDW phase, this expression differs from that in 
the paramagnetic state in two ways. First, the integrand 
contains velocities defined by ((SJ and, second, the in- 
tegration is performed over the modified Fermi surface, 
given by 



An analytic calculation is straightforward (see Ap- 
pendix 0, but here we estimate the integral iQ using 
geometric considerations, which will allow us to develop 
an intuitive picture. 

1. Crossing spheres 

First, we perform an estimate in the crossing case. 
When the gap opens, the original Fermi surfaces disap- 
pear in belts of width (ifci.2 and radius k± determined by 
the positioning of the spheres (see Figl^fA). For clarity, 
only Vi and dki are shown in the inset). The width of 
the disappearing belt can be evaluated using ((HJ. For an 
infinitesimal A, it implies that the original electron and 
hole Fermi surfaces arc destroyed in belts of width 

|vi X V2I ?l 

(see Appendix IXI for details). This formula is valid only 
in the crossing case, where vi x V2 7^ 0. Disappearance of 
the belts is accompanied by Fermi surface reconnection, 
i.e. new portions of the FS are created. As seen from 
Fig|2lA), the area of the new FS is also proportional to 
A. 

Before the opening of the gap, the belts were giving 
a contribution of approximately (w^ ^dki +v'^ 2'^fc2)27rfc^ 
to (Tqq. In this formula, the index a is the direction in 
which the conductivity is measured. For the spherical 
Fermi surfaces considered here, two distinct directions 
exist: one along the line connecting the centers of the 
FSi and FS2 spheres, and another one perpendicular to 
it. The conductivities along these directions are denoted 
cr|| and i7_L, respectively. For both of them, the velocity 
projections are zeroeth order in A, thus the contribution 
of the belts is proportional to their areas, i.e. first order 
in A. Fermi velocities on the reconnecting parts of the 
FS are given by (|SJ) and thus are also zeroeth order in A. 
Therefore the contribution of these parts is first order in 
A as well. 

Thus, when the belts disappear, and the reconnecting 
parts are created, the change of a is 

(5cr(lia) ~ A. (9) 

2. Touching spheres 

Next, we estimate the change of conductivity in the 
touching case. Here, (0J implies that the old Fermi sur- 
faces are destroyed in the circles centered at the touching 
point (see Figl^^B)). In a coordinate system with axes 
along Q and perpendicular to it, 

ei = vidku + 

" 2mi 

£2 = V2dk^\ + ^ (10) 



5 



That gives for both FSi,2 

dfc|| ^ dk']_ 

Using now © and (|10|l to obtain the radius dk± of the 
destroyed circle, we find dk±^ ^ VA- The area of the 
destroyed patches of the Fermi surfaces is A dk^ ~ A. 
The area of the reconnecting surface can be estimated 
as the area of a cylinder with a radius of order -s/A and 



height dki 



dk\ 



A. Thus the area of the reconnect- 



ing FS part is of order /S.^/"^. 

In contrast with the crossing case, the integrand of O 
in the touching case strongly depends on the direction a. 
For (T||, it is simply equal to [vl + i.e., zeroeth or- 
der in A, so the contribution of the destroyed part of the 
Fermi surface is {v\ + V2)A ~ A. By virtue of the argu- 
ment already used for the crossing case, the integrand for 
the reconnecting part is also zeroeth order in A. Since 
the area of the reconnecting Fermi surface is higher or- 
der then the one of the disappearing FS, its contribution 
can be neglected. Therefore, the change of the parallel 
conductivity is given by 



5(7 



(11) 



For (T_L, the integrand of Q is much smaller, because 
the projection of vp on the perpendicular direction van- 
ishes as one approaches the touching point. The char- 
acteristic value of the projection is vpidk/kp) ^ vp^/A, 
thus the contribution of the disappearing part of the FS 
to (7± is + i;|)AA ~ A^. For the same reason as in 
the calculation of a^, the value of the integrand on the 
reconnecting parts of the FS is of the same order \/A. 
Consequently, the contribution of these parts is of higher 
order (~ A^/^) and can be neglected. One obtains: 



da I 



A^ 



Formulae H9I11I12|I can be compared with those ob- 
tained by Elliott and Wedgwood^. For the touching case, 
our results are identical to theirs, while for the crossing 
case, Ref.!3 gives da\\ ^ A, da± ~ A^. The discrepancy 
with (jnj arises from the fact that Ref. H considers cross- 
ing of two spheres with equal Fermi momenta and Fermi 
velocities, while in the Rice model, generically, the sizes 
of the spheres differ. Coinciding Fermi momenta and 
Fermi velocities lead to the cancellation of tems linear in 
A, which cannot be captured by geometrical reasoning. 
This point is further explained in Appendix IXI 



B. Hall conductivity 



The change of this quantity upon onset of the SDW is 
found by calculating the difference between the integrals 
above in the SDW state (with the modified Fermi sur- 
face and Fermi velocities) and in the paramagnetic state. 
However, performing such a calculation for the linearized 
spectrum one finds Sa^^l^ = 0. This is due to the fact 
that although the curvature of the Fermi surface changes 
significantly upon reconnection of the two Fermi surface 
sheets, the two connecting arcs give equal contributions 
of opposite sign, which cancel each other exactly for a 
linear spectrum. Thus, one has to expand ei^2 up to 
quadratic terms in 5k even in the crossing case. 

Instead of performing the calculation based on H13|l . 
one can find the leading power in the dependence of 

( 

^"ap-y on A in a more intuitive way by using the ap- 
proach of Ong^"*. In this method, the Fermi surface is 
divided into slices perpendicular to the magnetic field 
B. The contribution of each slice to the Hall conduc- 
tivity is given by the area Ay inside the contour plot of 
v(yfc), k e FS. For B\\z, 



nil- 



where Ib = \/ ch/eB is the magnetic radius, and the total 
Hall conductivity is given by the integral 



'xy 



dkz A^(kz) 



27r 



nil 



(14) 



Variation of Ay{kz) at each slice kz upon onset of the 
SDW can be found from simple geometric considerations, 
which we describe below. 

In the Rice model, the Fermi surface has cylindrical 
symmetry around the direction of Q. Accordingly, there 
are two distinct elements of the tensor. We denote 



(12) them <t[^^ and crV"'', where parallel and perpendicular 



refers to the mutual orientation of Q and B. In the coor- 
dinate system where Q points along i, the conductivity 
tensor in a magnetic field has the form 



/ 



V 



^11 



Bz 'yrBy\ 



(15) 



( H) 

Below, we find the functional forms of a\"'{A) and 
(A). 



( H) ( H) 

The Hall conductivity cr-^ ' has the form a^^p 
'^ap-yB-y, where, for small magnetic fields—, 



dSp 



vp dkpdkg 



B^. (13) 



1. Crossing spheres 

Calculation of a^"' (B _L Q). Let the i-axis point 
along Q, and the z-axis point along B. The spherical 
electron and hole Fermi surfaces intersect along a cir- 
cle that lies in the plane kx = const. Following Ongii, 
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consider a cross-section of fc-space perpendicular to B 
(fcj = const cross-section). This cross-section cuts the 
electron and the hole Fermi surfaces along the circles and, 
depending on the value of kz, these circles may or may 
not cross. When they do cross, the contour plot of v is 
significantly modified due to opening of the gap. A typ- 
ical modification of the Fermi surface cross-section upon 
transition from the paramagnetic to the SDW state is 
shown in Fig. O as well as the corresponding modifica- 
tion of the V contour plot. Note that, in this plot, all k 
and V vectors are projections of the corresponding three- 
dimensional vectors on the = const plane. This means 
that the diameters of the circles in the contour plot of v 
are not vi and V2, but instead are given by the abso- 
lute value of the projections. However, to keep notation 
compact, we do not explicitly show this in the following 
discussion. 

Near the crossing points X in fc-space (Fig. O upper 
panel), the Fermi surface cuts reconnect. In the v-space, 
short segments disappear in the v contour plots, orig- 
inating from the electron and hole Fermi surfaces, and 
those contours reconnect as shown in the middle right 
panel of Fig. |3| The connecting segments are close to 
straight lines. Note that in formula (|14(l . the area inside 
the contour is positive or negative, depending on the con- 
tour orientation. Examining the reconnection pattern in 
Fig. 121 one concludes that the contribution da'^^^ under 
discussion can be estimated as the area difference of the 
two curved elongated triangles - ABC and CDE in the 
blow ups Fig.|ni;A,B). 

The triangle sizes depend on A. Let us first estimate 
the length of the the triangle bases dvl^2^ i-e. of the seg- 
ments that disappear from the original v-space contour 
plots. In the fc-space panel of Fig. O the length of the 
disappearing segments is given by JHJl. In the v-space, 
we first find the angles /3i.2 at which the gaps in FSi^2 
are seen from the centers of the respective Fermi sphere 
(see Fig. EI A)) 



|3^ 



dki 
kpi 



The lengths of the disappearing segments in the contour 
plots are 



dvi = Vifii = —^dki, 



vf 2A 
|vi X V2I nkpi 
vl 2A 



dvi 



Ivi X V2I hkFi 



Hence both triangles in Fig. 01 A,B) have areas propor- 
tional to A. The dominant contribution to the Hall con- 
ductivity would be (5(T^^^ ~ A, unless the triangles are of 
equal area, in which case only the higher order terms in 
A survive the subtraction. Here we show that this is not 
the case. Using the angles defined in Fig. |3IA,B), one 
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FIG. 3: Opening of the SDW gap modifies the Fermi surface 
(fc-space) and the contour plot of the velocity (v-space). In 
the velocity contour plots (before and after the SDW onset), 
the solid line corresponds to the electron and the dotted line 
to the hole surface. Panels (A) and (B) are blow-ups of the 
shaded rectangle on the left side of the v-space panel. They 
give the definitions of the angles Qi,2, /3i,2 and Q. Point O is 
the origin in the v-space. Points Gi,2 are the center points 
of the intervals dvi,2. Points A,B,D,E are defined as the 
points where deviations occur from the paramagnetic case. 
They are defined by the criteria Jsj. 



finds the ratio of the triangle areas 

dvi cos Qi cos ai fci?2f 1 



A2 



dv2COsa2 cosa2kpiV2 



(16) 



Due to the fact that ai, a2, and 6 are the angles in the 
triangle OG1G2 with the sides \OGi \ = vi and IOG2I = 
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V2, everything can be expressed through 9 



A2 



vi 

V2 + vi COS kpiV2 



(17) 



Thus, in the generic case, the two areas are different, and 
our conclusion about Sa^^^ ~ A holds. 

The next step is to integrate over k^. As mentioned 
earlier, some kz = const sections contain the crossing line 
of the electron and the hole Fermi surfaces, and some do 
not. The former sections contribute Sa^^^ ~ A, while 
the latter sections contribute a higher power of A, since 
these contours deform only slightly. In the case of cross- 
ing spheres, the number of contributing sections is pro- 
portional to fc_L, and thus is independent of A. As a 
result, the change of the whole integral H14|) is of the 
order of 



Calculation of a'l^^^ (B||Q). In this case, there are 
no contours that deform significantly, but a set of con- 
tours totally disappears. Each disappearing contour is of 
radius k±, and contributes 6a^^^ A ^ (/cx/fci?)'^ to the 
Hall conductivity. The main contribution to the integral 
over kx comes from the contours in the interval dk^ ~ A. 
Therefore, in the crossing case, and for BjlQ, one finds 



A. 



(18) 



Sa 



(H) 



A. 



(19) 



2. Touching spheres 



Calculation of CT^^' (B _L Q). The contour plot 
transformation differs from the case of crossing spheres 
in three ways. 

The first difference is that now the vanishing segments 
of the original contours dvi and dv2 are parallel to each 
other (in the notation of equation (|16|l . cosai — cosck2 = 
l). 

The second difference is that the area of each triangle is 
A ~ -s/A. This is because the length of the triangle bases 
dv is proportional to the length of the disappearing FS 
dk± and, as discussed in Sec. II VA^ and shown in Fig. |21 
dk± ~ VA. The areas of different sign still do not cancel 
each other in the generic case, so each section now gives 
a contribution dcr'^^' ~ \/A. 

The third difference is that now the number of sections 
with significant modification of the contours is also pro- 
portional to dkz = dk± ^ a/a. Applying the same logic 
as in the crossing case, we conclude that 
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( H') 

Calculation of cr[| ^ (B||Q). There are no contours 
that deform significantly, but a set of contours disap- 
pear altogether. Each vanishing contour had the radius 
dk± - VA and contributed da*-"^ ~ A ~ A. The mam 



(H) 



(20) 



contribution to the integral over k^ comes from the inter- 
val dkx ~ dk'^ ~ A (again, see discussion in Sec. II V A^ . 
As a result, in the case of touching spheres. 



A^. 



(21) 



C. Hall coefficient 



The resistivity tensor is the inverse of H15|) . Calculated 
up to the first order in B, it reads: 



1 




^ B, 
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B, 



B:r 



\ 



J 



Consequently, two distinct Hall coefficients can be de- 
fined. For B along Q, and current perpendicular to both 
of them, 



Sit) 



R 



and, for B perpendicular to Q, and current perpendicular 
to B, 



cr||(T_L 



Using equations [TTl [El UHl [H EH EH we get 



and 



6Rn 



(5i?|| 
SRj_ 



A' 
A 



(crossing case) 



(touching case) 



V. NUMERICAL RESULTS 

For the case of vanadium doped chromium, our ana- 
lytic results can be checked by doing a numerical calcula- 
tion based on the band theory Fermi surface of chromium. 
Such a calculation was presented in our earlier papeu^, 
where we emphasized the strong change in the Hall num- 
ber with doping, as observed experimentally^^. However, 
we also implied that the actual variation of the numerical 
results with doping near the critical point might be dif- 
ferent than that given by the analytic results based on a 
spherical Fermi surface. We did not, though, work close 
to the critical point because of convergence issues. This 
slow convergence can be traced to the problem of switch- 
ing from the antiferromagnetic e± basis to the paramag- 
netic £1.2 basis for the energy bands as the energy gap 
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goes to zero, which as we stated in the earher paper, 
would lead to a discontinuity in the Hall number for a 
perfectly flat Fermi surface. Since the actual Fermi sur- 
face always has some curvature, this problem can be com- 
pensated for by going to a more dense breakdown of the 
Brillouin zone when doing the momentum integration us- 
ing the linear tetrahedron method. 

We now describe how these calculations are performed. 
Chromium has three energy bands which cross the Fermi 
energy. One of them defines a hole octahedron around 
the H point as well as hole ellipsoids around the N points. 
The second defines an electron octahedron about the V 
point. The third defines electron balls centered around 
the T — H lines. (In this paper, we ignore the third band 
when computing the transport integrals, thus the elec- 
tronic structure becomes of the same two band form, as 
described in the previous sections.) The paramagnetic 
band structure is constructed from the local density ap- 
proximation using a linear muffin tin orbital method with 
muffin tin overlap corrections, and spin-orbit treated in a 
semi-relativistic fashion. After convergence, eigenvalues 
were generated at 506 k points in the irreducible (l/48th) 
wedge of the Brillouin zone. These were then fit by a 910 
function Fourier series (spline fit), which allows eigenval- 
ues to be generated for an arbitrary k. An SDW gap is 
then introduced as in the previous sections. Because of 
the convergence issues mentioned above, the calculation 
was carried out to even higher precision than in our previ- 
ous work. That is, the transport integrals are performed 
using the linear tetrahedron method with the momen- 
tum sums involving over 805 million tetrahedra in the 
paramagnetic zone. 

In Fig. 2{a) , we show results for the longitudinal ((T| | ) 
and transverse ((Tj_) conductivities (assuming a single Q 
domain). The variation is indeed linear in A for both 
cases, in agreement with the analytic results (crossing 
case). 

The A-dependence of ai^} is shown in Fig.^b). The 

value of aiyz hardly changes as the SDW gap A (equiv- 
alently, the magnetic order parameter) sets in. This is 
a consequence of the flatness of the octahedral surfaces 
along with the electron-hole compensation of the F and 
H sheets in the regions where they nest (the large ab- 
solute magnitude of ai^} is mostly coming from the N 
ellipsoids, which are essentially unaffected by the mag- 
netism). 

We can, though, test the analytic dependence derived 
earlier by looking at the individual contributions of the 
-I- and - bands, which we also plot in Fig. EJb). A clear 
linear A dependence is seen. Interestingly, the two linear 
terms have opposite sign with almost equal magnitude, 
which is a reflection of the electron-hole compensation in 
this basis. 

We have also tested our argument further by numeri- 
cally determining ai^} for the Rice model with unequal- 
sized spheres. Our numerical results indeed verify the 
analytic results presented earlier in this paper. More- 
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FIG. 4: (a) Calculation of a^^{axx) and (y±{<yyy) for chromium 
as a function of the energy gap, A (milli-Rydberg units). The 
chemical potential was set to correspond to four percent hole 
doping, with a Q vector of 27r/a(0.909, 0, 0). Empty symbols 
represent the result of a calculation in the single domain state. 
Filled symbols are the domain averaged result, (b) Calcula- 
tion of axyz for the domain averaged case. Although the small 
variation in each band (-1-,-) is linear in A, the signs are op- 
posite, leading to a nearly flat behavior of the total axyz with 
respect to A. 



over, we find that even in the paramagnetic case, axyz is 
small due to near electron-hole compensation, and also 
has a weak dependence on A for the same reason. As a 
consequence, axyz is primarily determined by the "reser- 
voir" , which stays essentially unchanged as the magnetic 
order sets in. Indeed, as we mentioned above, the domi- 
nant contribution to ai^} plotted in Fig.^b) comes from 
the N ellipsoids. 

As a consequence, we find the deviation of the Hall 
number to be linear in A as well, in agreement with the 
analytic results. Moreover, the near constancy of axyz 
implies that the Hall number should scale with the square 
of the longitudinal conductivity, a point discussed in our 
previous paper^. 

We have also checked the dependence of an on the 
variation q of the ordering wave vector Q, and found it 
to be of order q^, thus a higher order effect. We note 
that Q itself was chosen so as to yield the maximum 
in the susceptibility (that is, optimal nesting), the lat- 
ter calculated using the same linear tetrahedron method 
employed in the transport calculations. In addition, the 
occupation number change with A is small, meaning the 
chemical potential variation with A is also a higher order 
effect. Therefore, the linear variation of an and Rh with 
A we find both numerically and analytically appears to 
be a robust result. 



VI. CONCLUSIONS 

In conclusion, we have analytically estimated the lon- 
gitudinal and the Hall conductivities for a spin density 
wave metal. Our conclusions are two-fold. First, we find 
that, under generic circumstances, these transport coef- 
ficients depend linearly on the SDW energy gap in the 
small gap limit. These results have been confirmed by 
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numerical calculations using the actual Fermi surface of 
chromium. 

Second, we showed that the behavior of the Hall con- 
ductivity ai^I is special in the case of V-doped Cr in that 
the change Scr-xyl upon magnetic ordering is small. This 
is due to the fact that the nested regions of the hole and 
electron surfaces which are removed by magnetic order- 
ing give contributions which are approximately equal in 
magnitude but opposite in sign, as revealed by detailed 
numerical calculations based on the actual Fermi surface. 

The effect of this electron-hole compensation is partic- 
ularly obvious within the Rice model, where the critical 
mismatch of the electron and the hole surfaces is a small 
number of order A(0)/eir, A(0) being the value of the 
energy gap for perfect nesting, and ep being the Fermi 
energy. Hence, the total Hall conductivity of the two 
bands, and its change upon the gap opening, is reduced 
by a factor of order A(0)/eF compared with the contri- 
bution of each band. In both the case of the Rice model, 
and the true Fermi surface of chromium, the value of 
<Jxyz itself is largely set by the electron reservoir, i.e. by 
electrons that are unaffected by the magnetism. As a 
consequence, cri^i is nearly constant over a wide param- 
eter range across the quantum critical point. 

We now make a brief comparison with the 



experimenti 
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The latest pressure results in V- 



doped Cr— have provided access to the behavior in 
the immediate vicinity of the quantum critical point. 
The Hall coefficient was found to have a square root 
dependence on pressure: 5Rh ~ \/¥c~~P- The Neel 
temperature also has the same dependence and, since 
within mean field theory, A is proportional to T/v, 
we conclude that A ~ — P as well. Thus Rh 
indeed depends linearly on A, as we found theoretically. 
Experiment also finds that the Hall conductivity Oxyi is 
unchanged across the quantum critical point, and that 
Pxy scales with p^^, as we also expect. 

However, the longitudinal resistivity was found to vary 
as 6p ~ (pc — pY^'^.^ This is different from our theoretical 
prediction 5p ^ IS. ^ \pPc-~P and seems to be at vari- 
ance with the claim that pxy ~ p^^.. We believe this is 
connected with the difficulty of extracting critical expo- 
nents in the presence of sample inhomogeneity, though 
we cannot exclude the possibility of new physics beyond 
the mean- field picture. 

The overall success of the nested SDW picture for 
understanding the T = data in vanadium-doped 
chromium shows that the Hall effect is indeed a good 
measure of Fermi surface evolution across a quantum 
critical point. Our conclusion is not only important for 
simple metals such as chromium, but also serves as a 
benchmark for the understanding of the Hall effect in 
more complex systems such as heavy fermion metalsi^ 
and high Tc cupratesii. 

Finally, we should stress that our theory only applies at 
zero temperature, and does not address the strong tem- 
perature dependence of the Hall conductivity observed 



experimentally^^. We leave these interesting challenges 
for future theoretical studies. 



Acknowledgments 

We would like to thank Tom Rosenbaum, Gabe Aep- 
pli, Anke Husmann, and Minhyea Lee for communicat- 
ing their unpublished results to us. This work was sup- 
ported by the U. S. Dept. of Energy, Office of Science, 
under Contract No. W-31-109-ENG-38 (MRN, YB, RR) 
and by NSF Grant No. DMR-0090071, Robert A. Welch 
Foundation, and TCSAM (QS). 



APPENDIX A 



To simplify notation, we introduce 

i = (e i - £2)72 , 
Eo = Ve' + A2. 



Then 



ei + £2 



and, for the velocities in the modified bands, we find: 

V = i(vi-|-V2)±^(vi-V2), (Al) 

where — Vi-^. Since the integral in Q is taken over 
the modified Fermi surface, we first have to evaluate ex- 
pression IjAip at the new FS. We start with the electron 
FS (equation |(SJl with the "+" sign) for which (|A1|) is 
equivalent to 



(A2) 



^ ^ ( 1 + I; I , 

, Eo 

where the factors u'^,v'^ are mathematically identical to 
the BCS coherence factors. At the Fermi surface, one 
can express all quantities through ei using Q). Keeping 
in mind that for the electron FS ei < 0, we calculate 



En = 




— (A3) 
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Now one can express the u^^v^ factors at the FS: 



Thus, for the electron FS, 
A2 



rvi 



rV2, 



(A4) 



(A5) 



AV + 2A^ej(vi.V2)+e>i 
(e2 + A2)2 



Note, however, that m this equation vi_2 must be evalu- 
ated on the FS, so if the velocities vary substantially, the 
usage of ljA5|l is not straightforward. 

For the hole FS there is a sign change in ljA2p . Also 
now ei > and there is an additional sign change in the 
formula l|A3|l for Eq. As a result, the formula for v comes 
out to be the same. Thus Eq. ljA5|l is valid for both parts 
of the FS. 

The expression for the conductivity |(7J| is written in 
cartesian coordinates {k^, ky, kz). However, it is much 
easier to evaulate it in the (ei, £2, (j)) coordinates, (j) being 
the angle of rotation around the Q axis, because of the 
simplicity of Eq. Q in these coordinates. To do that, 
we need to express the element dSp of the Fermi surface 
area in the (ei, £2, (f>) coordinates. 

Fig.l^A) shows a section of the electron and hole Fermi 
spheres in a plane cutting through their centers, i.e. a 
(j> = const section. All (f> = const sections are identical 
in our model. In the plane, we introduce cartesian coor- 
dinates {k\\,k±) with fc|| being along Q. Using a linear 
approximation for the energies £1,62 near the point X, 
where FSi and FS2 cross, one has 



ei 

£2 



vi • (?idk), 
V2 • {hdk), 



(A6) 



with vi,2 taken at point X . It is convenient to introduce 
two reciprocal vectors Ui,2, such that 



Explicitly, 



Ul 



U2 



Uj = S,j. 



V2 X [Vi X V2] 

(Vi X V2)2 
Vi X [V2 X Vi] 



(Vi X V2)2 



In terms of Ui 2, one has 



A2 , dei 



dk — Uidci + U2C?e2 = (ui 



(A7) 



(A8) 



where the second expression is valid only at the Fermi 
surface. Squaring equation IjASI) . we obtain 

^^2 ^ vi + 2(AVef)(vi.V2) + (AVef)vf / de,^ ' 

(vi X V2)2 V n 

(A9) 



Using IjASp and IIA9|l together, we find 

dk _ ej + A'^ jdei/n) 

\vf\ e? |vi X V2I ' 

dSp e? + A2 (dei/n) 

1 r = 2 — I \kj_d(j), 

\vf\ i-i |vi X V2I 

where k\_ is the radius of the circle over which the two 
original FS cross. If the integrand has no </) dependence, 
one can write 



dSp 

M 



27rfc_L 
|vi X V2I 



-|-oo 



dei t\ + A2 



(AlO) 



Here integration over the (—00, —0) interval accounts 
for the electron FS contribution, and integration over 
(+0, -l-cx)) accounts for the hole FS contribution. 

We now continue with the calculation of the conduc- 
tivity tensor iTa/3. With our choice of coordinates, it is 
diagonal, and its two unequal components are g\\ and o ^. 
To use integral O, note that 



2 / 2 I 2 \2 

Wll = [U Will -^V V2\\) 



^±1 



(u^WlJ. -I- W^W2_l)^ 



(All) 



Start with g\\. Using ljA4|l in IjAlip we can now write 
down the expression for the integrand of Q. Up to a 
constant factor (?t j (AifiK) one gets: 



2-Kk±_ 

Vi X V2 I 

dei t\ + A2 



27rfc 



Vi X V2 



U WJii + 2u W t'i||U2|| + V 1'2|| f ~ 

(A12) 



" / X^fTA^{f^^ll+2AV|P2||+e?.2^||}. 

Notice that this integral does not converge near = 0. 
This only means that our approximation l|A6|l does not 
allow us to calculate the total value of cr, since it is accu- 
rate only in the vicinity of the crossing point X . However 
we can calculate the change of the conductivity, ba . To 
do that, we need to subtract the result corresponding to 
the paramagnetic case from (IA12|I : the sum of the inte- 
grals over the original FSi,2 (i.e. along the dashed lines 
on Fig. HA)). 

Up to the same constant factor e^T j {Aic'K) that we 
chose not to show explicitly in ljA12|l , the contribution of 
FS2 in the paramagnetic state is 



,(2) 



27rfc 



Vl X V2 



de\ 



^2| 



where deiX)2/(?i|vi x V2I) = dk-i is exactly the length ele- 
ment of FS2. By the same token, the FSi paramagnetic 
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contribution is 



APPENDIX B 



.(1) 



We get 



5a\\ 



27rfc 



Vi X V2 



<h2 

h ''ill 



2TTk± f dci 



Vi X V2 



27rfci 



|vi X V2I 



^7rfc^(2wi||U2|| - u^ii - 1^211) 



^ 2 2 
~T^l\\ ~ ""211 
'1 



Vi X V2 



+°° &i 

-00 ^ e? + A2 



27rfc_L(wi|| - ^211)^ A 



|Vi X V2I 



(A13) 



Note the factor with the difference of Fermi velocity pro- 
jections in this formula. For the special case of iden- 
tical Fermi surfaces, it may vanish for some directions, 
and then the change of longitudinal conductivity will be 
quadratic in A. This indeed happens in the case consid- 
ered in Ref. and was discussed in the body of the paper 
at the end of Sec. lIV A2l In the Rice model with different 
Fermi surfaces, such a cancellation does not happen. 

For the other components of the conductivity, we 
merely need to use the corresponding velocity projec- 
tions in An additional complication here is the 0- 
dependence in (|A11|) . but it only gives a factor of 1/2: 



Vi X V2 



7:k^{vii_ - V2i_)^ A 



Vi X V2 



The conclusion for the crossing case is that 



6a 



(l|:-L) 



A 



(A14) 



(A15) 



In the case of a linear SDW, each sheet of the FS in- 
teracts with the other one shifted both by -l-Q and — Q. 
For an incommensurate Q, this leads to an infinite sec- 
ular matrix which has to be solved. On the other hand, 
for Fermi surface restricted integrals, we can formulate a 
good approximation which avoids this difficulty. For the 
Rice model in the incommensurate case, for instance, it 
never happens that the copies of the hole Fermi surface 
shifted by -l-Q and — Q are both close to the same point 
of the electron Fermi surface. This property consider- 
ably simplifies the calculation, because it allows us to 
consider the hole FS shifted only by -|-Q. Indeed, the 
electron Fermi surface can be divided into three parts: 
(I) the one which disappears due to interaction with the 
hole Fermi surface shifted by -l-Q, (II) the one which dis- 
appears due to the interaction with the other hole surface 
shifted by — Q, and (III) the one which survives. Any ki- 
netic coefficient will be given by part III. At the same 
time, when only one of the shifted surfaces is considered, 
it is given by a sum of parts II and III for -l-Q, and parts 
I and III for — Q. Since the paramagnetic value is a sum 
of the contributions from I, II, and III, it follows that the 
linear SDW leads to coefficients which are equal to the 
sum of the values for the -|-Q and the — Q cases minus the 
paramagnetic value. Therefore, the change of the kinetic 
coefficients is twice the result obtained in the calculation 
with the hole surface shifted by only -fQ, since the -l-Q 
and — Q contributions are equal by inversion symmetry. 



A similar argument can be applied to the ffat Fermi 
surface model of Shibatani et aM^, as well as to the actual 
Fermi surface of Cr. In Section V, the above relation was 
considered when doing the numerical calculations as well. 



1 E. Fawcett, Rev. Mod. Phys. 60, 209 (1988). 

^ E. Fawcett, H. L. Alberts, V. Yu. Galkin, D. R. Noakes 
and J. V. Yakhmi, Rev. Mod. Phys. 66, 25 (1994). 

^ M. Lee, A. Husmann, T. Rosenbaum and G. Aeppli, un- 
published. 

* P. Coleman, C. Pepin, Q. Si and R. Ramazashvili, J. Phys.: 

Condens. Matter 13, R723 (2001). 
^ Q. Si, S. Rabello, K. Ingersent, and J. L. Smith, Nature 

413, 804 (2001); Phys. Rev. B68, 115103 (2003). 
^ J. A. Hertz, Phys. Rev. B 14, 1165 (1976). 
^ A. J. Millis, Phys. Rev. B 48, 7183 (1993). 
® R. J. Elliott and F. A. Wedgwood, Proc. Phys. Soc. London 

81, 846 (1963). 

^ M. R. Norman, Q. Si, Ya. B. Bazaliy and R. Ramazashvili, 

Phys. Rev. Lett. 90, 116601 (2003). 
" S. Chakravarty, C. Nayak, S. Tewari and X. Yang, Phys. 

Rev. Lett. 89, 277003 (2002). 
^ J. Ziman, "Electrons and Phonons", Ch.l2 (Oxford Uni- 



versity Press, London, 1960). 

A. Yeh, Y.-A. Soh, J. Brooke, G. Aeppli, T. F. Rosenbaum 
and S. M. Hayden, Nature 419, 459 (2002). 
T. M. Rice, Phys. Rev. B 2, 3619 (1970). 
N. P. Ong, Phys. Rev. B 43, 193 (1991). 
Here we have assumed that the change of r due to the 
magnetic ordering can be ignored. This is expected to be 
valid for intermediate values of the impurity scattering po- 
tential. It is also experimentally found^ to be the case for 
the pressure-tuned transition in Cri-^^Vi: 
S. Paschen, T. Luehmann, C. Langhammer, O. Trovarelli, 
S. Wirth, C. Geibel, F. Steglich, Acta. Phys. Pol. B 34, 
359 (2003). 

F. F. Balakirev, J. B. Betts, A. Migliori, S. Ono, Y. Ando, 
and G. B. Boebinger, Nature 424, 912 (2003). 
A. Shibatani, K. Motizuki and T. Nagamiya, Phys. Rev. 
117, 984 (1969). 



